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Abstract 

The lifetimes of the B° and B~ mesons are measured using a sample of about 
four million hadronic Z decays collected from 1991 to 1995 with the Aleph 
detector at LEP. The data sample has been recently reprocessed, achieving a 
substantial improvement in the tracking performance. Semileptonic decays of 
B° and B~ mesons are partially reconstructed by identifying events containing 
a lepton with an associated D* + or D° meson. The proper time of the B 
meson is estimated from the measured decay length and the momentum of 
the D-lepton system. A fit to the proper time of 1880 D* + £~ and 2856 D°£~ 
candidates yields the following results: 

r B o = 1.518 ± 0.053 ± 0.034 ps, 
r B - = 1.648 ± 0.049 ± 0.035 ps, 
t b -/t bo = 1.085 ±0.059 ±0.018. 
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1 Introduction 



Measurements of the individual b hadron lifetimes represent an important test of the 
present knowledge of nonspectator effects in the b hadron decay dynamics, such as Pauli 
interference, W exchange and weak annihilation. Based on the heavy quark expansion 
formalism, the difference between the lifetimes of the b baryons and mesons is predicted 
to depend on I /ml and higher order terms, whereas meson- meson differences depend only 
on 1/mf and higher order terms [|IJ. The predicted hierarchy is r\ b < t b o ~ t Bs < t b - . 
Differences are expected to be at the level of a few percent, which sets the scale of the 
experimental precision required. 

This paper reports an improved measurement of the B° and B~ lifetimes with the 
Aleph detector at LEP, using approximately four million hadronic decays of the Z, 
collected in the period 1991-1995. This data sample was recently reprocessed, achieving 
higher efficiency and better resolution in the track reconstruction, which is highly 
beneficial for this analysis. 

Semileptonic decays of B° and B~ mesons are partially reconstructed by identifying 
events containing a lepton (electron or muon) with an associated D° or D* + meson. The 
resulting D°-lepton (D°£~) and D* + -lepton (D* + £~) event samples consist mostly of B~ 
and B° decays, respectively (charge conjugate modes are implied throughout this paper). 

Previous measurements of the B° and B~ lifetimes are reported in |2|. |3|. 

2 The ALEPH detector 

A detailed description of the Aleph detector can be found elsewhere j|, ||. A high 
resolution vertex detector (VDET) consisting of two layers of silicon with double-sided 
readout provides measurements in the r<p and z directions at average radii of 6.5 cm 
and 11.3 cm, with 12 /mi precision at normal incidence. The VDET provides full 
azimuthal coverage, and polar angle coverage to |cos#| < 0.85 for the inner layer only and 
Icos^l < 0.69 for both layers. Outside the VDET, particles traverse the inner tracking 
chamber (ITC) and the time projection chamber (TPC). The ITC is a cylindrical drift 
chamber with eight axial wire layers at radii of 16 to 26 cm. The TPC measures up to 21 
space points per track at radii between 40 and 171 cm, and also provides a measurement 
of the specific ionization energy loss (dE/dx) of each charged track. These three detectors 
form the tracking system, which is immersed in a 1.5 T axial magnetic field provided by a 
superconducting solenoid. The combined tracking system yields a momentum resolution 
transverse to the beam axis of <j(pt)/pt = 6 x 10~ 4 pr © 0.005 (pt in GeV/c). The 
resolution of the three-dimensional impact parameter for tracks having two VDET hits 
can be parametrized as a = 25 /mi + 95 [im/p (p in GeV/c) 0. 

The electromagnetic calorimeter (ECAL) is a lead/ wire-chamber sandwich operated 
in proportional mode. The calorimeter is read out in projective towers that subtend 
typically 0.9° x 0.9° in solid angle, segmented in three longitudinal sections. The hadron 
calorimeter (HCAL) uses the iron return yoke as absorber. Hadronic showers are sampled 
by 23 planes of streamer tubes, with analogue projective tower and digital hit pattern 
readout. The HCAL is used in combination with two double layers of muon chambers 
outside the magnet for muon identification. 
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Recently the LEP1 data were reprocessed using improved reconstruction algorithms. 
The features that are particularly relevant for the enhancement of the charmed meson 
reconstruction efficiency are the following. A new VDET pattern recognition algorithm 
allows groups of several nearby tracks which may share common hits to be analysed 
together, to find the hit assignments that minimize the overall x 2 f° r the event. The 
improvement in the hit association efficiency is more than 2% (from 89.2% to 91.0% in 
r(f) and from 85.6% to 88.2% in z). Information on the drift time from the TPC wires 
is combined with that obtained from the pads to reduce the error in the z coordinate 
by a factor of two. A 30% improvement in the rip coordinate resolution is achieved for 
low momentum tracks by correcting the pad coordinates for ionisation fluctuations along 
the tracks as measured by the wires. The particle identification (dE/dx) is improved by 
combining pulse height data from the TPC pads with that of the wires. The improvements 
in the charmed meson reconstruction efficiencies with respect to the previous analysis 
range from 10 to 30%, depending on the decay channel. 

3 Event selection 

The D* + £~ and D°£~ event samples consist of an identified lepton (e or fi) associated 
with a fully reconstructed D* + or D° candidate. The selection of muons and electrons 
is described in detail in [0. For this analysis, lepton candidates are required to have a 
momentum of at least 2.0 GeV/c for electrons and 2.5 GeV/c for muons. 

The D* + and D° candidates are reconstructed from charged tracks and 7r°'s that 
form an angle of less than 45° with the lepton candidate. These charged tracks are also 
required to intersect a cylinder of radius 2 cm and half-length 4 cm centered on the 
nominal interaction point, to have at least 4 hits in the TPC, a polar angle 9 such that 
|cos#| < 0.95 and a transverse momentum greater than 200 MeV/c. 

Photons and 7r°'s are identified in the ECAL. The four- momenta of 7r° candidates are 
computed by adding the photon momenta when the 77 invariant mass is consistent with 
the 7T° mass. The energy of the tt° is then recomputed using the kinematical constraint 
of the 7T° mass ||. The energy resolution achieved is about 6.5%, almost independent of 
the energy. For reconstructing the D° candidates, only 7r°'s with momenta greater than 
2 GeV/c are used. 

For all charged kaons used to reconstruct the D° it is required that the dE/dx be 
within 3cr of that expected from a kaon, except for the D* + — > D°tt + , D° — > i^~7r + decay 
channel. The kaon has the same sign as the lepton coming from the B semileptonic decay, 
therefore this charge correlation is required. 

Kg candidates are reconstructed from pairs of oppositely charged tracks. The two 
tracks are required to be inconsistent with originating from the interaction point and the 
Kg candidate is rejected if the measured mass is more than 2cr (±10 MeV/c 2 ) from the 
nominal K§ mass. 

Tracks coming from the D° decays are required to form a common vertex with a x 2 
confidence level greater than 1%. In the D° decay channels that do not contain a 7r° 
in the final state, if more than one combination satisfies the selection criteria, the one 
with the smallest x 2 °f the D vertex fit is selected. For the decay modes where the 
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multiple candidates originate from different re combinations, for a given detected lepton, 
the candidate with m^o^ — m D o closest to the nominal value is chosen in case of D* 
selection, otherwise the highest momentum D° candidate is chosen. 

In order to reject the background from charm, to improve the signal to background 
ratio and to ensure well-measured decay lengths, additional selection criteria are applied 
to all the subsamples. The invariant mass of the D^*'£ system is required to be greater 
than 3 GeV/c 2 , where D^*' can be D* + or D°. This cut significantly reduces the charm 
background while keeping ~ 85% of the signal. To exploit the high precision of the silicon 
vertex detector, the lepton track and at least two tracks from the D° decay are required 
to have at least one VDET hit in both the r<p and z projections. 

3.1 The D*-lepton selection 

D* + candidates are identified via the decay D* + — > D°tt + , followed by D° — > K~tt + , 
~7T + 7T 7T + , "7r + 7T° or 7r + 7r . The difference in mass between 

the D* + and D° candidates must lie within 1.5 MeV/c 2 (approximately two standard 
deviations of the experimental resolution) of the nominal value of 145.4 MeV/c 2 . 

In the D° — > K~tt + channel, the D° momentum p^o is required to be greater than 
5 GeV/c. 

For the D° — > K~7t + 7t~tt + channel pr>o must be greater than 8 GeV/c, and at least 
two of the D° decay tracks must have p > 1 GeV/c. 

In the case of the D° — > K^7i + 7T° channel, the momentum of the reconstructed D° is 
required to be greater than 10 GeV/c, and the two charged tracks in the decay are required 
to have p > 0.5 GeV/c. Futhermore the decay kinematics are required to be consistent 
with one of the three resonant decays: D° — > K~p + , D° — > K*~tt + , D° — > K*°7c°. For 
each decay hypothesis the mass of the resonant particle and the helicity angle 9jj are 
calculated^] If the mass is consistent with the nominal resonance mass value, within its 
natural width, and if | cos##| > 0.4, the candidate is considered to be consistent with the 
resonant decay hypothesis. For these resonant decays a cos 2 9h distribution is expected. 

Finally, in the D° — > K®it + tt~ channel, the momentum of the reconstructed D° is 
required to be greater than 4 GeV/c. The same technique as in the D° — > K~7t + tt° 
channel is used to tag the K*~ resonance. The sign of the resonance is used to distinguish 
between D° and D°. 

3.2 The I^-lepton selection 

The D°£~ sample consists of events with a lepton and a D° candidate, where the D° is 
not the decay product of a reconstructed D* + . The D° candidates are identified using 
the same decay modes as for the D* + £~samp\e. For this sample, the background is larger 
because the D* + -D° mass difference criterion is not applicable. 

For all the decay modes selected in the D°£~ sample, a search for the additional pion 
is performed to reject D° candidates coming from D* + — > D°tt + . If a pion candidate is 

1 The helicity angle 9h is defined as the angle between the scalar particle and one of the decay products 
of the vector particle, calculated in the rest frame of the vector particle. 
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found yielding a D* + -D° mass difference within 6 MeV/c 2 of the nominal value, the D°£~ 
candidate is rejected. The efficiency for reconstructing the additional pion and rejecting 
D 0, s coming from D* + decays is found to be 86% from the Monte Carlo simulation. 

In the D° — > K~n + channel, the reconstructed D 0, s are required to have 
p D o > 8 GeV/c, p K > 2 GeV/c and p„ > 1.5 GeV/c. In the D° — > i^7r + 7r~7r + channel 
the .D momentum is required to exceed 12 GeV/c while the kaon and the three pions 
must have momenta greater than 2 and 1 GeV/c, respectively. 

The cuts on the kaon and pion momenta are tightened for the D° — ► K~7t + 'k 
channel, to p% > 3 GeV/c and > 2 GeV/c, while the reconstructed D° must have 
p D o > 12 GeV/c. The same cuts on the three resonances are used as in the .D* + ^~event 
sample. 

Finally, for the D° — > i^g7r + 7r~ channel a cut is applied of at least 1.5 GeV/c on 
the kaon momentum and 1.0 GeV/c on the two pion momenta. The momentum of the 
reconstructed D° candidates is required to be greater than 10 GeV/c. The same technique 
is used to tag the K*~ resonance as in the D* + £~ event sample. 

3.3 B meson reconstruction 

The B decay vertex position is estimated by vertexing the reconstructed D° track with 
the lepton. Events are rejected if the B vertex fit gives a x 2 probability less than 1%. 

The D° candidate mass spectra for the four subsamples in the D* + £ event selection 
are shown in Fig. |], and in the D°£~ event selection in Fig. |[ The fitted curves consist 
of a Gaussian function for the signal plus a linear background. For the D° — > K~7t + tt° 
in the D°£~ sample a Gaussian tail is used to describe the background; masses below 
1.65 GeV/c 2 are excluded to avoid the broad enhancement due to the missing tt° in 
the D° — > K~7t + 7t°tt° decay. The fitted D° width and the fitted number of signal and 
background events within a window of ±2<r around the fitted mass are shown in Table [I| 
for the different samples. 

Events reconstructed within two standard deviations of the fitted D° mass are selected 
for the lifetime analysis, resulting in 1880 D* + £~ and 2856 D°£~ candidates. The decay 
length is calculated for these events by reconstructing the primary and B decay vertices in 
three dimensions. The primary vertex reconstruction algorithm M applied to simulated 
bb events yields an average resolution of 50 /zm x 10 /im x 60 (ttrn (horizontal x vertical 
x beam direction). 

The distance between the primary and B decay vertices is projected onto the direction 
defined by the momentum of the D^£ system. The uncertainty on the flight direction 
due to the missing neutrino induces a negligible error on the decay length. The resolution 
on the B decay length is on average 250 /im, compared with an average B decay length 
of about 2.5 mm. 

For the D° — ► K~7t + 7t° channel, the tt° momentum is included when extrapolating the 
neutral D° track backwards to form the B vertex. In the case of D* + £~ events, the soft 
pion from the D* + decay does not improve the resolution on the B decay length and is 
therefore not used in the reconstruction of the B vertex. 

Because the selected decays contain an undetected neutrino, the B momentum is not 
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Figure 1: The invariant mass of D° candidates for the four subsamples in the D* + £ event 
selection: a) D° — > K~ir + , b) D° — > K~ir + ir~ir + , c) ~7T + 7r° (notice the different mass 

scale), d) D° — ► K^tt^tt + . The superimposed curves are the results of the fit described in the 
text. 

known precisely and is reconstructed using an energy flow technique as described in ||. 
A further correction is applied by evaluating for Monte Carlo events the k distribution, 
defined as 

097) 

reco 1 \ 

Because this distribution depends on the selection criteria applied, separate k distributions 
are calculated for each subsample. The momentum resolution obtained is between 8 and 
10%, depending on the decay channel. 
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Figure 2: The invariant mass of D° candidates for the four subsamples in the D°£ event 
selection: a) D° — > K~tt + , b) D° —* K~7r + ir~7r + , c) 'tt + tt° (notice the different mass 

scale), d) D° — > i^g7r + 7r~. The superimposed curves are the results of the fit described in the 
text. 

4 Lifetime measurement 

An unbinned likelihood fit for the lifetimes is performed. For each event the probability 
of observing a proper time t given the lifetime r is calculated: the probability density 
function F(t, a t , r) is obtained by convoluting an exponential distribution with the 
properly normalised k distribution and with a Gaussian function which takes into account 
the resolution on the decay length. 

Both the D* + £~ and D°£~ samples contain a mixture of B° and B~ decays and the 
B~ j B° mixture in the samples depends on the ratio of the lifetimes, as discussed in Section 
[4.2| . Therefore, to measure the B° and B~ lifetimes a simultaneous maximum likelihood 
fit is performed to all the events. The likelihood function contains three components for 
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Table 1: Fitted D° width (a), number of D° candidates and fraction of background events 
falling within a mass window of ±2<r. The uncertainties are statistical only. 



Subsample 


Width (MeV/c 2 ) 


Candidate events 


Background fraction 




D° - 


-> K-71+ 


8.6 ±0.3 


651 


0.066 ± 0.004 




D° - 




6.6 ±0.3 


670 


0.096 ±0.004 




D° - 


-> K-TT + TT 


30.3 ± 1.6 


394 


0.127 ±0.008 




D° - 


-> K^TT+TT' 


8.3 ±0.7 


165 


0.061 ±0.006 


D°£r 


D° - 




8.7 ±0.3 


1312 


0.133 ±0.006 




D° - 




7.2 ±0.5 


664 


0.232 ±0.012 




D° - 




28.4 ± 1.9 


563 


0.258 ±0.012 




D°- 


-> Kgvr+Tr- 


8.6 ±0.8 


317 


0.139 ±0.009 



each sample and is written as 



x 



f[ [f-(TB-/T B o)F(t t , a t , t b -) + ^{r B -/r m )F{t h a 4 , r B o) + /4 g Kk g &) 
i=i 

II [f-(T B -/TB°)F(U, <n, r B -) + f^r B -/r B0 )F(t u a h r B o) + f^F^U)} . (2) 



i=i 



The coefficients /* and /q are the fractions of the D* + £~ sample arising from B~ and B° 
decays, respectively. Similarly, f°_ and /q are the fractions of B~ and B° decays in the 
D°£~ sample. The coefficients /£ kg and /° kg are the background fractions of the samples, 
while the functions F£ k (t) and F® k (t) are their normalised proper time distributions. 

4.1 Backgrounds 

Background contamination arises from the following sources: 

(1) combinatorial background, i.e. candidates with a fake 

(2) the process B — > DjD^X, followed by Dj — > £~X, giving rise to a real and 
a real lepton; 

(3) a real meson accompanied by a fake or nonprompt lepton, from Z — > bb or 
Z — > cc events. 

Source (1) is the dominant background. Its contribution is determined from a fit to the 
D° mass distributions, and its magnitude is given in Table |l| for the various subsamples. 
The proper time distribution for this source is determined from the data by selecting 
events from the high mass sideband of the D° peak. The same selection criteria described 
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in Section |3] are applied to the background samples, except that the requirement on 
the D* + -D° mass difference in the case of the D* + £~ events is removed to increase the 
statistics. A function consisting of a Gaussian plus positive and negative exponential tails 
is used to describe these data. 

The contribution from source (2) is calculated from the measured branching ratios 



for this process |L0| and a Monte Carlo simulation to determine the detection efficiency; 
it accounts for a contamination which is about 2-3% of the sample, depending on the 
channel. 

The background from source (3) is estimated from the measured hadron-lepton 
misidentification probabilities and the measured inclusive D° and D* + rates. An 
independent estimate is obtained using wrong-sign (D* + £ + or D°£ + ) events, and is 
found to be consistent. This background source contributes between 2 and 5% of the 
sample, depending on the channel. To characterize the proper time distribution of this 
background, it is further subdivided into three distinct components: fake leptons coming 
from the primary vertex in cc and bb events, and fake leptons coming from a decaying b 
hadron. 

The proper time distribution for sources (2) and (3) are determined from simulated 
events. 



4.2 Sample compositions 

Both the D* + £~ and D°£~ samples contain a mixture of B° and B~ decays. In order 
to estimate the cross contamination, the individual semileptonic branching ratios of the 
B° and B~ must be determined. The evaluation follows the same procedure as given 
in the appendix of 0. An important input in this evaluation is the knowledge of the 
branching ratios for the decay modes B — > D**^^, where D** represents any of the 
p— wave resonances as well as the nonresonant D^ir states. The most recent ALEPH 



and DELPHI results JTT|, [12| for both the resonant and the nonresonant components are 
used in the calculation, leading to a significant reduction in the resulting uncertainty 
compared to the previous analysis of 0. The B° and B~ content of the two samples are 
calculated using as input the measured values of the branching ratios given in Table 
The sample composition is then calculated by considering the B° and B~ decay channels 
that contribute to the D* + £~ and D°£~ samples ||, |14| , taking into account the probability 
of 0.147 ± 0.015 that a D* + £~ event is mistakenly reconstructed SIS db D°£~ event. 

As a consequence of this procedure, the coefficients /*, /q, f°_ and /q appearing in 
the likelihood function (Eq. depend on the lifetime ratio. For equal lifetimes, about 
85% of the B decays in the D* + £~ sample are attributed to B°, while about 80% of the 
D°£~ sample B decays come from B~ . 

4.3 Fit results 

The fit to the proper time distributions of the D* + £~ and D°£~ events is performed to 
determine the two free parameters rgo and r B - . The values obtained are 

t B o = 1.518 ± 0.053 ps, 



8 



ALEPH 




Figure 3: Proper time distributions with the result of the fit overlaid for the two samples. 
The plots show the background contributions to the samples, together with the B° and B~ 
components. 

r B _ = 1.648 ± 0.049 ps, 

where the errors are statistical only. The statistical correlation coefficient is —0.35. The 
ratio of the lifetimes is 

t b -/tbo = 1.085 ±0.059, 

taking into account the correlation. 

The proper time distributions for the two samples are shown in Fig. |3], with the results 
of the fit overlaid. 
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4.4 Systematic uncertainties 



The sources of systematic uncertainty are discussed in the following, and the estimated 
errors are summarized in Table |3|. 

The uncertainty in the B momentum reconstruction is dominated by the uncertainty 
in the k distribution. The effects that can modify the k function with a consequence on 
the extracted lifetimes are studied. The first is due to the D^'tt content of semileptonic B 
decays, which affects the k function due to the smaller phase space available and the softer 
momentum spectrum of the reconstructed D^i system. The fractions of D^it decays 
are varied in the simulation within experimental errors. Other effects studied are the 
dependence of the k distribution on the analysis cuts and the b fragmentation function. 
These two effects combined give a relative uncertainty of about 2% in the B momentum 
determination. Uncertainties in the k distribution propagate in the same way to tb- and 
Tgo, and therefore have a small effect in the ratio. 

Uncertainties in the background fractions and proper time distributions are considered. 
The combinatorial background, estimated in the data from fits to the D° mass spectra, is 
reported in Table |l] for each individual channel. The systematic uncertainty due to this 
source is estimated by varying the combinatorial background in the fit for each channel 
within its statistical errors in turn and taking the sum in quadrature. The contamination 
from physics background is evaluated using the simulation, by varying the fraction within 
the estimated uncertainty and repeating the fits. The uncertainty on the total background 
level is calculated by combining the two above sources. 

The parameters describing the background proper time distributions are varied within 
their uncertainties. Background proper time distributions are parametrized using different 
methods, to check for possible systematic bias. Different background samples are selected 
by varying the sideband regions, adding events from the lower sideband, which are 
excluded in the lifetime determination, or using events with wrong-sign correlations. Some 
cuts in the selection are varied to check the stability of the parametrizations. The shapes 
extracted from real data are compared with those extracted from Monte Carlo events. The 
resulting differences in the fitted lifetimes are used to estimate the systematic uncertainty. 

For the remaining systematic errors, which are small compared to the statistical ones, 
the correlation between the lifetimes is not propagated into the ratio. 

The systematic uncertainty due to the sample compositions is determined by varying 
the branching fractions of Table ^| within ±lcr from the central values. In addition the 
uncertainty due to the assumption of isospin conservation has been estimated, allowing 
for a variation of 20% relative to the exact symmetry. 

As explained in detail in 0, a small fraction of four-body decays where an extra pion 
is produced in the decay of the B meson contributes to the sample composition, and is 
characterized by a different selection efficiency with respect to the real signal. The relative 
efficiencies 

-±=- , ; = 0.75 ± 0.10 (for B°) 0.64 ± 0.10 (for B~) 

e(B — > D^'iu) 

therefore enter into the calculation and this uncertainty (coming from the limited Monte 
Carlo statistics) is propagated to the measured lifetimes. 
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Table 2: Branching ratios used as input values in the calculation of the sample composition. 



Decay 


B.R. 


Reference 


B° -> D+e-D 


0.0460 ±0.0027 
0.0200 ±0.0025 
0.1045 ±0.0021 


i 
1 


13 

1 

13 


! 


B° -»■ D°n+£~iJ + 
5° -> D*°7r + £-u 


0.0159 ±0.0036 


11 


, 12| 


Br -> D*+7r-£-v 


0.0121 ±0.0018 


|11 


, 12| 


_> D+7r-£-p 


0.0124 ±0.0048 
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B~ -> D + ti~£~v + 
B~ -> D* + tt-£-u 


0.0157 ±0.0031 


[ 
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Table 3: Sources of systematic error on the fitted lifetimes. 



Source of error 


Contribution to systematic error 


t B o (ps) 


tb- (ps) 




B momentum reconstruction 


±0.025 


±0.026 


±0.009 


Background treatment 


±0.020 


±0.020 


±0.010 


Sample compositions 


±0.003 


±0.003 


±0.004 


D(*'ir£~is relative efficiency 


±0.006 


±0.006 


±0.006 


Decay length resolution 


±0.008 


±0.008 


±0.008 


Total 


±0.034 


±0.035 


±0.018 



The parameters of the decay length resolution function are varied within their errors 
to estimate the resulting uncertainty [14|] . 



5 Conclusions 



The lifetimes of the charged and neutral B mesons have been measured with the full 
statistics collected by the Aleph detector at and around the Z peak energy. The data 
sample was recently reprocessed, achieving improved tracking performance. A maximum 
likelihood fit to the proper time distributions of 1880 D* + £~ and 2856 D°£~ candidates 
yields the following results for the B° and B~ lifetimes and their ratio: 

t bo = 1.518 ± 0.053 ± 0.034 ps, 
r B - = 1.648 ± 0.049 ± 0.035 ps, 
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tb- /tb° 



= 1.085 ± 0.059 ± 0.018 , 



where the first error is statistical and the second is systematic. 

These results supersede the corresponding ones of 0. Averaging with the results of the 
other methods presented in p| the combined values are: 
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